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Abstract 

We present a theory of magnetization reversal due to thermal fluctuations in thin submicron-scale 
rings composed of soft magnetic materials. The magnetization in such geometries is more stable 
against reversal than that in thin needles and other geometries, where sharp ends or edges can 
initiate nucleation of a reversed state. The 2D ring geometry also allows us to evaluate the effects 
of nonlocal magnetostatic forces. We find a 'phase transition', which should be experimentally 
observable, between an Arrhenius and a non-Arrhenius activation regime as magnetic field is varied 
in a ring of fixed size. 
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I. INTRODUCTION 



The dynamics of magnetization reversal in submicron-sized, single-domain particles and 
thin films has attracted much attention given its importance in information storage and other 
magnetoelectronic applications. The problem can be approached by stochastic methods: in 
the classical regime (typically at temperatures above ~ 1°K) the magnetization dynamics 
is governed by the Landau-Lifschitz-Gilbert equation 111 perturbed by weak thermal noise. 

nn 

The classical Neel-Brown theory (2J, |3|] of thermally induced reversal assumed a spatially 
uniform magnetization and uniaxial anisotropy. Experimental confirmation of this theory 
has been provided for certain simple single-domain systems (15-30 nm Ni, Co, and Dy 
nanoparticles) Q. 

There nevertheless remain fundamental open questions, especially when there is spatial 
variation of the magnetization density 0, B, Isj] . While in small particles that are spherical 
or nearly so, as in j^, the Neel-Brown theory appears to work reasonably well, it appears to 
break down for elongated particles, thin films, and other geometries, which exhibit far lower 
coercivities than predicted 

Braun [5J] made an initial step by studying the effects of spatial variation of magnetization 
density on magnetic reversal in an infinitely long cylindrical magnet. However, Aharoni 
pointed out that the energy functional employed neglected important nonlocal magneto- 
static energy contributions, invalidating the result. Further, for submicron- scale magnets 
with large aspect ratio, finite system effects are likely to play an important role; for exam- 
ple, simulations [f], ^| indicate that magnetization reversal in cylindrical-shaped particles 



proceeds via propagation and coalescence o 
ends. Both of these issues are addressed in 



magnetic 'end caps', nucleated at the cylinder 
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Here we consider a geometry that avoids these difficulties: an effectively two-dimensional 



annulus. Such systems have recently received increasing attention |13l . Il4j . They are typi- 
cally constructed of soft magnetic materials (quality factor Q ~ O(10~ 2 )), such as Fe, fee 
Co, or permalloy, have radii of order (10 2 — 10 3 ) nanometers and thicknesses of order 10 nm 
or less. 

Our interest in these systems is twofold. The first is technological: because the mag- 
netic bending length is much smaller than the typical system size, there are two oppositely 
polarized stable states, each with magnetization vector pointing everywhere along the cir- 
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cumferential direction; they are degenerate in the absence of an external magnetic field. But 
a current running along the z-direction through the center, with z the direction normal to 
the annulus plane, generates a circumferential magnetic field breaking the degeneracy. By 
switching the direction of the current, the relative stability of the two states is switched. 
(A slightly different method, but with similar wire dimensions and current magnitudes, 
was used in The utility of such a system as an information storage device depends 

on the magnetization state being relatively long-lived against thermal fluctuations, even at 
relatively high temperatures. Unlike the cylindrical particle, the micromagnetic ring has 
no ends where nucleation is easily initated, making its magnetic state more stable against 
thermally induced reversal. 

The second is physical: by developing a theory for thermally induced reversal that can 
be solved analytically, we are able to extract a number of interesting qualitative features 
that would be more difficult to uncover numerically and which should apply also to more 
complicated situations. While several important quantitative features require a numerical 
treatment, we show below that our most important qualitative findings are robust (see in 
particular the Discussion section). 

One of our central predictions is that a type of phase transition occurs in the thermally 
induced reversal rate, and more importantly, that it can be realized experimentally. The 
possibility of such a transition in classical stochastic field theories (which we show below 
includes the physical problem of interest here) was first noted in |l5l |. as system size was 
varied in a symmetric Ginzburg-Landau double-well 4 potential. It was further shown in 
to apply more generally to asymmetric systems as well. In the present case, the transition 
depends on two parameters: the system size and the strength of the applied magnetic field. 
Although the former cannot be continuously varied, the latter can, facilitating experimental 
tests of the predicted transition. In particular, we will show that as magnetic field varies for 
a ring of fixed size, there should be a transition from a regime where activation is Arrhenius 
to one where it is non- Arrhenius. 

A preliminary account of this work has appeared in 
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II. THE MODEL 



We consider an annulus of thickness t, inner radius R\ and outer radius i?2- We confine 
our attention to rings satisfying t <C AR <C R, where AR = R 2 — R± and R = + R^jl. 
A current run through the center leads to an applied field H e at R in the circumferential 
direction 9; the small variation ~ 0(AR/R) of field strength with radius can be ignored. 
As will be discussed below, magnetostatic forces produce strong anisotropics, forcing the 
magnetization vector to lie in the plane and preferentially oriented parallel to the inner and 
outer circumferences. We may therefore consider magnetization configurations that vary 
only along the ^-direction. 

Suppose now that the system is initially in its metastable state; i.e., with magnetization 
vector M = —MqO. We are interested in determining the mean rate for thermal fluctuations 
to reverse the magnetization to its stable direction. We consider temperatures above 1°^, 
where classical thermal activation can be expected to apply. The magnetization dynamics 
are then governed by the Landau-Lifshitz-Gilbert equation 

<9 t M = - 7 [M x H cff ] + (a/M )[M x d t M] , (1) 

where M Q is the (fixed) magnitude of M, a the damping constant, and 7 > the gyro- 
magnetic ratio. The effective field H e ff = — 5E/5~M. is the variational derivative of the total 
energy E, which (with free space permeability /io = 1) is 
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18]: 



£[M(x)] = X 2 [ d 3 x\VM\ 2 + - [ d 3 x\VU\ 2 

- / d 3 xH e -M, (2) 



Jn 

where Q is the region occupied by the ferromagnet, A is the exchange length, |VM| 2 = 
(VM,) 2 + iyM y f + (VM 2 ) 2 , and U (defined over all space) satisfies V ■ {VU + M) = 0. 
The first term on the RHS of Eq. (j2j) is the bending energy, the second the magnetostatic 
energy, and the last the Zeeman energy. Crystalline anisotropy terms are neglected, given 
their negligibly small contribution; they can be easily included but will at most result in a 
small modification of the much larger shape anisotropies, to be discussed below. 
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III. ENERGY SCALING AND THE MAGNETOSTATIC TERM 



The presence of the nonlocal magnetostatic term complicates analysis. However, the 
quasi-2D nature of the problem allows a significant simplification, as shown by Kohn and 
Slastikov in an asymptotic scaling analysis that applies when the aspect ratio k = t/R 
and the normalized exchange length I = X/R are both small, and I 2 ~ k \ log A; | . These con- 
straints restrict the range of ring geometries to which our analysis applies. Before discussing 
the KS result, we recast the energy in dimensionless form, letting X = x/R and similarly 
for all other lengths, and let h = H e /(2M / 2 ). Then, integrating along the direction normal 
to the plane, we have for the bending plus Zeeman energy contribution 

E b + E z 



M 2 i? 3 



kl d 2 X 



(V x m) 2 -2h-m , (3) 



where m = M/Mo is the normalized magnetization and u represents the 2D surface with 
boundary duo. 

Before analyzing these terms further, we examine the magnetostatic energy contribution. 



The analysis of 
terms: 



12 1 showed that this asymptotically separates into local bulk and surface 



JpL = hjj 2 Xm 2 z + (l/4n)k 2 \\o g k\ l^xim-rf 

+ U? j d 2 X|V-m| 2 ff _ 1/2 , (4) 

where sx is dimensionless arc length along the boundary, and the final integral is the squared 
H^ 1 / 2 Sobolev norm of V-m j^j]. With current technology, the orders of magnitude k ~ 10~ 2 
and I 2 ~ k\ \ogk\ ~ 10~ 2 — 10 _1 are just attainable. Then the first term of Eq. (J3J is larger 
than the others by roughly two orders of magnitude, forcing m z = 0, and therefore in this 
topology we can ignore fluctuations of m out of the plane (we will discuss this further in 
Sec. EH}. 

The second term, like the first, is a (local) magnetostatic surface (or shape anisotropy) 
term. The third term represents a nonlocal magnetostatic bulk energy. When nonzero, this 
term will be roughly an order of magnitude smaller than the others (see Sec. IV BJ) . and so 
to a first approximation 2l( it can be dropped. This will result, for some values of ring size 
and external field, in an error of up to 10% in the computation of the action, so we can hope 
at best for reasonably good quantitative predictions of the logarithm of the escape rate. As 



FIG. 1: Ferromagnetic annulus viewed from above, showing coordinates used in text. 

noted in the Discussion section, however, the important qualitative features uncovered by 
our analysis should remain unaffected. 

We use a locally varying coordinate system where the angle <p(s r ) measures the deviation of 
the local magnetization vector from the local applied field direction; i.e., = indicates that 
the local magnetization is parallel to the local field, <fi = ir indicates that it is antiparallel, 
and so on. The parameter s' = R9 is the arc length along the circumference. The geometry 
and variables used are displayed in Fig. [T] 

The normalized magnetization vector can therefore be written, in cylindrical coordinates, 
as m = (m r ,mg,m z ) = (sin0, cos0, 0). After integrating out the radial coordinate the 
bending plus Zeeman energy becomes 



The m r -^ — term is a 'winding number' of with respect to the local direction; 

it gives zero in all configurations considered here, but would give a nonzero contribution 
for uniform magnetizations, e.g., m = x. For fixed Mq, it does not contribute to the 
magnetization equations of motion given by Eq. ((TJ. 

Finally, subtracting out constant terms and the first derivative term (which gives zero 
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contribution), noting that the boundary integral occurs over both inner and outer radii, and 
rescaling lengths again gives 



£ 



1/2 



ds 



{— ) 2 + sin 2 <f> — 2h cos< 

OS 



(6) 



where £ = E/E = E/[2M$R 2 AR y /c'kl% s = 6y/c, £ = 2ity/c, h = |H e |/(2M / 2 c), and 
c = (1/2tt)(A;| log k\/l 2 )(R/AR). In deriving © we used the fact that log(R 2 /Ri) = AR/R+ 
O (AR/R) 2 . The error is negligible for the geometries considered here; for example, with 
the ring parameters used in Fig. [7[ AR/R = 1/5 and log(-R 2 /-Ri) = -20067. The parameter 
c (0 < c < oo) depends on the ring size and material properties; it represents the ratio of the 
anisotropy energy scale to the bending energy scale, and determines the width of a Bloch 
wall. 



A. Energetics and Topology 

The scaling results of the previous section are useful insofar as they provide results on how 
different contributions to the energy scale in the thin- film limit. Their effective application in 
a physical situation must also take into account the geometry, and in our case, the topology 
of the ferromagnetic particle under study. Consideration of both of these aspects provides 
a guide for considering what types of magnetization configurations might be relevant in 
different thin-film geometries. 

As one example, we consider the flat disk topology studied by Shinjo et al. [22] (cor- 
responding to our geometry with Ri = 0). They studied magnetization configurations in 
permalloy disks of thickness 50 nm and diameters ranging from 0.3 to 1 micron. Inter- 
estingly, they observed vortex structures, particularly in the larger diameter samples. The 
surface term energy in Eq. (£Q) (which leads to the shape anistotropy) is minimized by requir- 
ing the magnetization vector to remain tangential to the surface (i.e., in the i^-direction). 
But given the topology of the samples used here, this forces the interior magnetization to 
do one of two things: either the magnetization magnitude goes to zero at the center, or 
the magnitude stays mostly constant but then the out-of-plane magnetization component 
m z ^ in some interior region. Either choice costs energy, but (when considered over the 
same region) the first costs more than the second. 

Given that there must be an out-of-plane magnetization component in the disk topology 
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studied by Shinjo et al. |22|, the KS analysis can determine the approximate lengthscale 
over which this occurs. Comparing the first term in Eq. (J3J) with the bending term (j3J), 
we estimate that their respective energies are of the same order when the the region where 
m z 7^ is roughly of the order of an exchange length, i.e., 10-20 nm. This appears to 
be exactly what is observed. (Shinjo et al. [22] don't provide an estimate for the width 
of this region, noting only that they observe a contrast 'spot' at the center of each disk 
that corresponds to out-of-plane magnetization; see Fig. 2 of their paper.) Note that if the 
magnetization were to be out-of-plane in a region much larger than this, Eq. (0J) predicts a 
prohibitively large energy cost. (Simulation results consistent with these conclusions appear 
in |23|.) 

In our ring topology, however, an 'outer vortex' configuration (i.e., magnetization circum- 
ferential at the outer boundary) does not require an out-of-plane magnetization anywhere. 
We can therefore ignore configurations with m z ^ 0, any of which are likely to have energies 
larger than the configurations considered here. 

IV. TRANSITION IN ACTIVATION BEHAVIOR 

The reversal rate T due to thermal fluctuations at temperature T is given by the Kramers 
formula T ~ Fo exp(— AW/ksT), where the activation barrier AW ^> ksT is simply the 
energy difference between the (meta)stable and 'saddle' states. The latter is the state of 
lowest energy with a single negative eigenvalue (and corresponding unstable direction) of the 
linearized zero-noise dynamics. Equiyalently, it is the configuration of highest energy along 



the system's optimal escape path 24]. The rate prefactor r is determined by fluctuations 
about this optimal path, and its evaluation will be presented in Sec. IV CI 

Stable, unstable, and saddle states are all time-independent solutions of the LLG equa- 
tions. For fixed Mo, Eq. (0) and the variational equation H c fj = — 5E/5M. yield a nonlinear 
differential equation that must be satisfied by any such time-independent solution: 

d 2 (p/ds 2 = sin ^ cos + /isin0. (7) 

There are three 'constant' (i.e., (ft is independent of 6; these remain nonuniform configurations 
because m varies with position) solutions for < h < 1: the stable state = (m = 6); the 
metastable state (ft = n (m = — 9), and a pair of degenerate unstable states (f> = cos -1 (— h), 
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which constitute the saddle for a range of (£, h). The = 0, 7r solutions are degenerate when 
h = 0, and the <ft = ir solution becomes unstable at h — 1. We therefore confine ourselves to 
fields in the range < h < 1. 

We have also found a nonconstant ('instanton') solution j^j of Eq. (|7J), which we will see 
is the saddle for the remaining range of (£, h). It is 



{s,s ,m) 



2 cot" 1 



s - so, \ sn(7£|m) 
dn — - — m 1 



(8) 



5 ) cn(JZ\m) J 

where dn(-|m), sn(-|m), and cn(-|m) are the Jacobi elliptic functions with parameter m, 
< m < 1 [26(; So is an arbitrary constant arising from the rotational symmetry of the 
problem; and 1Z and 5 are given by 



sn 2 (ft|m) = 1/m- h/2 - (l/2m)^m 2 h 2 + 4(1 - m) (9) 
5 2 = =- . 10 

2 - (m + ^2/12 + 4(1-771)) 

The period of the dn function equals 2K(m), the complete elliptic integral of the first 
kind (2^. Accordingly, imposition of the periodic boundary condition yields a relation 
between t and m: 

£ = 2K(m)5. (11) 

In the limit m — *■ 1, corresponding to £ — » 00 at fixed /i, Eqs. fl%j)- (fTTj) reduce to Braun's 
solution P|. In the limit m — > 0, dn(x|0) = 1, and the instanton solution reduces to the 
constant state = cos _1 (— h). At m = 0, the critical length and field are related by 



4 = tt5 c = 2tt/ VI - K . (12) 

The solution (J%|)- (fTTj) corresponds to a pair of domain walls of width 0(5). At fixed h, 
the constant configuration is the saddle for i < i c and the instanton is the saddle for £ > £ c . 
This can be understood as follows: at fixed field, the bending energy becomes sufficiently 
large at small £ so that the constant state becomes energetically preferred. (There is a 
second transition at even smaller £, where the bending energy becomes so large that the 
magnetization lies along a single Euclidean direction everywhere; we do not consider such 
small length scales here.) Conversely, at fixed £ the constant configuration is the saddle for 
h > h c , and the nonconstant for h < h c . (However, when £ < 2tt, the constant configuration 
is the saddle for all h.) Here the Zeeman term dominates at sufficiently large field, preferring 
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FIG. 2: The phase boundary between the two activation regimes in the (£, /i)-plane. In the shaded 
region the instanton state is the saddle configuration; in the unshaded region, the constant state. 



a constant configuration, while at smaller field the shape anisotropy energy dominates, 
preferring the instanton configuration. The 'phase boundary' (Eq. ([12))) is the m = line in 
the (£, h) plane, and is shown in Fig. |21 We now compute the reversal rate in both regimes, 
and examine how it is affected by the transition in the saddle state. 



V. THE REVERSAL RATE 



We turn now to a computation of the magnetization reversal rate T due to thermal 
fluctuations at temperature T. In equilibrium, it is given, as noted in Sec. II VI by the 

n 

Kramers formula T ~ To exp(— AIF/fcsT) [24]. We first compute the activation barrier AW 
for each saddle configuration. 



A. Activation Energy 

As noted earlier, the exponential dependence of the magnetization reversal rate on tem- 
perature is given by AW, the energy difference between the saddle (<p u ) and metastable 
(0 S ) states (the notation arises from the properties that the saddle is unstable along the 
longitudinal escape direction, while the metastable state is locally stable in all directions). 
With the latter given by 4> s = ir, this is 

AW/E = £[0 u ]-£[0 s = tt] (13) 
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12.5 15 



FIG. 3: Activation energy AW for fixed h = 0.3 as I varies. The dot indicates the transition from 
constant to instanton saddle configuration. 



C/2 



ds 



;^) 2 + sin 2 (0 u )-4/ i cos 2 (0 u /2) 

OS 



When the constant state <p = cos _1 (— h) is the saddle configuration, it easily seen that 
AW = (l — h) 2 £/2. When the nonconstant, or instanton, state is the saddle, the integral (JT3)) 
must be computed numerically. However, it can be analytically computed in the m — > 
(£ — > t^{h)) limit, where one finds AW(m — > 0) — > (1 — h) 2 £/2. So the energy (and its 
first derivative, which can also be computed) is continuous at £ c {h). Of course, the second 
derivative is discontinous there. Fig. El shows the activation energy as a function of ring 
circumference for fixed field. 

The activation energy grows linearly with I when the transition state is constant; it 
becomes almost flat above l c at fixed h, because the width of the domain walls remains 
essentially constant (cf. Fig 4 of ^|). In Fig. 0] we show the activation energy dependence 
on h at fixed £, which is more relevant to experiment. 



B. Bulk Magnet ost at ic Energy Contribution 

We can now go back and check whether the contribution of the bulk magnetostatic term is 
small compared to the bending energy. This requires an evaluation, or at least an estimate, 
of the H- 1 ' 2 Sobolev norm of the divergence of the reduced (i.e., in-plane) magnetization. 
To do this, we need to introduce some additional notation. The L 2 (u) norm of a quantity 
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h) for L=7 




0.2 



FIG. 4: Activation energy AW for fixed £ = 7 as h varies. The dot indicates the transition from 
instanton to constant saddle configuration. 



(say, the gradient of the reduced magnetization), is 



I Vml 



l 2 h 



/ d 2 x (Vm) : 



1/2 



(14) 



so that the dimensionless bending energy is simply kl 2 \\ VmH^^)- From here on we shall 
abbreviate L 2 (u) to L 2 for ease of notation. 

Formally, the H^ 1 / 2 Sobolev norm of the magnetization divergence is given by 



\V-m\\ H -y* = ||(V)- 1/2 (V-m)|| L2 . 



(15) 



Its meaning becomes clearer through the use of Fourier transforms. Define the Fourier 
transform /(£) of /(x) in the usual way: /(x) = J dP£f(£) e^' x . Then 

2 



IV • ml 



d 2 H 



£irhi + 6^2 



1/2 



2t {€imi + &m 2 ) 2 



lei 



It now follows in a straightforward fashion that 



IIV7 II / II II / || 1 1 1/2 ||V7 || 1/2 

|| V • m | |zr-i/2 < ||m||^i/2 < ||m|| L2 ||Vm|| L2 



< ||Vm 



1 1/2 



(16) 



(17) 



where the last inequality follows because ||Vm|| L 2 = 1. Therefore, the bulk magnetostatic 
term is dominated by the bending energy. 
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As noted earlier, the relevant scaling regime for the approach presented here corresponds 
to £ 2 ~ k\ logk\ ~ 10~ 2 — lO^ 1 . For the constant saddle configuration the maximal bulk 
magnetostatic energy arises when U = ir/2; this is also the maximum value of <fi for the sad- 
dle, corresponding to h — 0. For this configuration, and with ring parameters used in Fig. 
an upper bound for the magnetostatic bulk energy, computed using the inequalities (|17|). is 
roughly an order of magnitude smaller than the bending energy. As h increases from 0, and 
correspondingly <p u — > ir, the magnetostatic bulk term decreases to zero. 

For the nonconstant, or instanton, saddle, the minimum value of i is 2tt. At this length- 
scale, an upper bound for the ratio of magnetostatic bulk energy to bending energy varies 
roughly from 0.05 to 0.1 as h varies; smaller numbers are found as lengthscale increases, 
justifying the neglect of this term. 

Qualitatively, the instanton configuration has nonzero divergence only over a region of 
0(5), which remains smaller than 0(1) except close to m — 1 (£ — > oo) and h = 1. The 
instanton configuration contributes to the bending energy, however, over the entire ring. It 
is therefore not surprising that, in the appropriate scaling region, the instanton's magneto- 
static bulk energy is relatively small compared to its bending energy. This is in contrast 
to instanton configurations in the cylinder P, 1(3]; there, while the region contributing to 
a bulk divergence is 0(1), the same region supplies the entire contribution to the bending 
energy as well, and so the magnetostatic contribution cannot be neglected there. 

C. Rate Prefactor 

The leading-order rate asymptotics are determined by the activation barrier AW; the 
subdominant asymptotics appear as the rate prefactor r . Because the magnitude of T 
is controlled by the extent of fluctuations about the optimal escape path, the prefactor is 
considerably more difficult to calculate than AW. Although the reversal rate is only linearly 
dependent on the prefactor, as opposed to its exponential dependence on the activation 
barrier, the rate can still be significantly affected by T , especially in the vicinity of a 
transition in saddle configurations (cf. Fig. |SJ). Moreover, an understanding of the prefactor 
is needed to study other quantities of physical interest, such as exit location distributions 
The prefactor computation procedure is summarized in 



24, 



28| (see also |2J,|3(|). Con- 



sider a small perturbation rj about the metastable state, so that sufficiently close to it 
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<f) = (ps + rj. Then to leading order the time dependence of fluctuations about the metastable 
state is given by rj — —A s i], where A s is the linearized zero-noise dynamics at 4> s . Similarly 



A u is the linearized zero-noise dynamics around <p u . Then 

r = 



24. 




det A s 



(18) 



det A u 

where A Uj o is the single negative eigenvalue of A u . Its corresponding eigenvector is the 
direction along which the optimal escape trajectory approaches the transition state. In 
general, the determinants in the numerator and denominator of Eq. ()18|) can separately 
diverge: they are typically products of an infinite number of eigenvalues with magnitude 
greater than one. However, their ratio, which can be interpreted as the limit of a product 
of individual eigenvalue quotients, is finite. 



1. Constant saddle 

When t < £ c (h), or equivalently h > h c (£), the saddle is the pair of constant configurations 
<p = cos _1 (— h), the prefactor can be determined by direct computation of eigenvalues of the 
stable and unstable states Q|. Linearizing around the stable state gives 

f] = -A s r] = -{-d 2 /dx 2 + l-h)r], (19) 

and similarly 

7] = -A u r] = -(-d 2 /dx 2 + h 2 - 1)77 (20) 
about the transition state. The spectrum of eigenvalues corresponding to A s is 

47T 2 r) 2 

K = -^ + ±~h n = 0,±l,±2... (21) 

and the eigenvalues corresponding to A u are 

4-ir 2 r) 2 

K = ^^ + h 2 -l ra = 0,±l,±2... . (22) 

This simple linear stability analysis justifies the claims that (f) s is a stable state and 4> u 
a saddle. Over the interval [0, £ c ) all eigenvalues of A s are positive, while all but one of A u 
are. Its single negative eigenvalue Aq = —(1 — h 2 ) depends on h but is independent of I. 

Putting everything together, we find 

r - = r -i (1 ~ h 2 ) sinh(yT^M/2) 
T ° vr sin(vT^7^/2) ' 1 ' 
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FIG. 5: The prefactor Tq (in units of Tq vs. h when I = 7 on the 'constant saddle' side of the 
transition. The prefactor on the 'instanton saddle' side of the transition acquires an additional 
temperature dependence, as discussed in the text. 

where Tq 1 = ajEo/MoV(l + a 2 ), with V the ring volume. The rate includes a factor of 2 
because the system can escape over either of the saddles, which are rotationally equivalent 
with respect to (f) s = it. 

The prefactor Tq diverges at £ c (h), or conversely h c (£), as expected (cf. Fig. EJ); in this 
limit, Tq ~ const x (£ c — £)~ l as i — > £~ at fixed h, or as (h — h^ 1 as h — *■ at fixed i. The 
prefactor in this region for fixed I as h varies is plotted in Fig. The divergence arises from 
the vanishing of the eigenvalue of a pair of degenerate eigenfunctions at the critical point. 
This indicates the appearance of a pair of soft modes, resulting in a transverse instability of 
the optimal escape trajectory as it approaches the saddle. The meaning and interpretation of 
the divergence is discussed in detail in Q]. Near (but not at) the critical point the prefactor 
formulae hold, but in a vanishing range of T as £ c is approached. Exactly at £ c the prefactor 
is finite but non- Arrhenius (with a different exponent than that in Eq. (|30|l). Inclusion 
of higher-order fluctuations [3l| about the saddle can be used to compute the prefactor at 
criticality, and will be addressed elsewhere. We return to the prefactor divergence in Sec. I VII 

The independence of Tq with respect to temperature leads to the well-known exponential 
temperature dependence of the overall reversal rate. By analogy with chemical kinetics, this 
exponential falloff of the rate is often called 'Arrhenius behavior'. 
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2. Nonconstant (instanton) saddle 



Computation of the determinant quotient in Eq. ()18|) is less straightforward when the 
transition state is nonconstant, i.e., when £ > £ c (h) or equivalently h < h c (£) . An additional 
complication follows from the translational degeneracy (energy invariance with respect to 
choice of so) of the nonconstant state. This implies a soft collective mode in the linearized 
dynamical operator A u of Eq. (jl8J) . resulting in a zero eigenvalue for all h < h c (£) (not to be 
confused with the vanishing of the lowest stable eigenvalue of the saddles exactly at h c (£)). 

To proceed, we use the McKane-Tarlie regularization procedure |32j, which allows the 
evaluation of det' A u , the functional determinant of the operator A u with the zero eigenvalue 
removed. We refer the reader to for details, but sketch the main features here. 

With periodic boundary conditions, it is formally the case that 

det 7 A u = 772(3 + t, s ; m) - i] 2 (s, s ; m) 
(VilVx) r} 1 {s,s Q ;m)detH(s,So;m) 

where r]i(s, s ;m) and r] 2 (s, s ;m) are two linearly independent solutions of A u r]i = 0, 

i = 1,2, (771I771) = f_i/ 2 ds yf(s,0;m) is the square of the norm of the zero mode and 

detH(s, So;m) = r) 2 (s, s ; m )Vi( s i s o'i m ) ~ Vi( s i s o'i m ) r l2(s, s ; m) is the Wronskian; here a 

dot denotes a derivative with respect to s. The expression f!24j) is meaningful only as part 

of a determinant quotient, as noted above. 

The functions rji and r] 2 can be found by differentiating the instanton solution (jHJ) with 

respect to s and m, respectively; i.e., r]i(s, s ; m) = d(fi(s, s ; m)/ds and rj 2 (s, So;m) = 

d(f)(s, s ;m)/dm. This yields 

i]i(s, So;m) = — sn(TZ\m)cii(TZ\m) 



sn(^ 2fi |m)cn( £: ^ fl |m) 
x ; i ; —, —, 9, — 25 



and 



cn 2 (^|m) + sn 2 (^|m)dn^(^o|m 



r] 2 {s,s ;m) = — - — — - x 

cn 2 (7c|m) + sm(7c|m)dn (^-^|m) 

m(s - s ) d5 f s-So s - s 

— : — sn( l-i mjcnf l-i mjsnf — - — mjcnf — - — \m) 

o 2 dm 

sn( — - — |mjcn( — - — \m)L{ — - — \m) 
000 



^sn(7?.|m)cn(7^|m) 



2(1 -m) 

/S — So, . ,S — So, v, S — So- 2 / S — S 1 Nj /S — So, v 

— sn( — - — mjcnf — - — mHl — m){ — - — ) — sn ( — - — mjdnf — - — \m) 
00 000 
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+dn(7?.|m)- — dn( — r~~\ m ) 
dm o 



} (26) 



where E(-|m) is the incomplete elliptic integral of the second kind 
Inserting these solutions into Eq. (ffifj) yields 

det' A u 



S 3 [(2m/6)(d6/dm)K(m) - K(m) + E(m)/(1 - m)) 



(27) 



4m 2 sn (1Z \ m) cn (7Z \ m) dn (7Z \ m) dlZ/dm 

Using a similar procedure, we find the corresponding numerator for the determinant ratio 
in (HHJ) to be 

det A s = 4 sinh 2 (6y/l - hK(m)) , (28) 

consistent with the numerator of Eq. ()23j) . obtained through direct computation of the 
eigenvalue spectrum. (Recall, though, that it is only the ratio of the determinants that is 
sensible.) This becomes clearer by noting that the expressions in Eqs. ()27|) and (|28J) are 
well-behaved for all finite i > £ c (m > 0). While both expressions separately diverge as 
m — > 1, it is easily checked that the divergences cancel. 

As already noted, the rotational symmetry of the instanton state (corresponding to the 
arbitrariness of the constant sq) corresponds to a 'soft mode', resulting in appearance of a 
zero eigenvalue A Uj i = of the operator A u . The corresponding eigenfunction is clearly rji 
given by ()25|) . The appearance of a zero mode corresponds to the zero rotational energy 
of the instanton solution: the center of the domain wall pair can appear anywhere in the 
ring. This is in contrast to the situation in the finite cylinder J^, where the instanton is 
'pinned'. The general procedure for including the correction required due to removal of the 
zero eigenvalue is described by Schulman J^iJ. The correction in our problem is an additional 
factor of 2£a/ /i^ksT ', which vanishes as m — > (and thereby removes the divergence 
of the pref actor as m — > + ). 

Finally, we need to compute the eigenvalue A^o corresponding to the unstable direction. 
With the substitution z = (s — s )/5, the eigenvalue equation A u r] = Xi] becomes 



(A + 1 - h)cn 4 {Tl\m) + 2(A + 3)sn 2 (ft|m)cn 2 (ft|m)dn 2 (z| 


m) + (A - 1 - h)sn\n\m)dn\z 


m) 


{cn 2 {TZ m) + sn 2 (7?. m)dn 2 (2 


\m)) 2 



The lowest eigenvalue corresponds to a nodeless solution for r\. By continuity it must tend 
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lambda (h) for L=7 



-0 . 65 
-0 . 675 
-0.7 
-0 . 725 
-0 . 75 
-0.775 
-0.8 




0.1 0.2 0.3 0.4 



FIG. 6: Lowest eigenvalue A Uj o as a function of h for I = 7. 

towards — 1 + h 2 as m — ■> + . We have solved Eq. (J2Hj) numerically for I = 7; the result 
appears in Fig. |H1 The weak dependence on h (and also £) is typical. 

Finally, we put all of the above results together to find the formula for the prefactor per 
unit length: 

T T+/i = \X (£, h)\m(k B T)- 1/2 sinh UVT^KKim) 

x 2mK(m)^ + ^[E(m) - (1 - m)K(m)ll "'^ (30) 
dm 1 — m 

where E(m) is the complete elliptic function of the second kind 2(|. As noted above, Xq(£, h) 
is weakly dependent on h and £, and is 0(1) everywhere. 

The most important qualitative feature to be noted from Eq. (fHUJ) is that the zero eigen- 
value arising from the uniform translation mode leads to non-Arrhenius behavior — i.e., a 
T-dependent prefactor — everywhere on the low-field side of the transition. 

Finally, we note that the eigenfunction rji given by (|25|) has a single pair of nodes. Because 
nodes arise in pairs, there must then be only a single (nodeless) solution of lower eigenvalue 
than rji. But rji has zero eigenvalue, proving that the solution (jHJ has a single unstable 
eigenmode, and is therefore a proper saddle. 

The above results allow one to find the overall reversal rate in any part of the (£, h) phase 
plane. Results for a permalloy ring with given dimensions are shown in Fig. [7| 

Among commonly used soft ferromagnetic materials, permalloy has the largest magnetic 
exchange length. The discussion of scaling in Sec. 11111 suggests that the effects of nonlocal 
magnetostatic terms are minimized with larger exchange lengths. Where else might one find 
magnetic materials with large exchange lengths? Such materials would require both low 
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FIG. 7: Total switching rate (in units of s 1 ) vs. (3 = 1/kgT (in units of °K 1 ), at fields of 
(a) 60 mT (instanton saddle) and (b) 72 mT (constant saddle). Parameters used are k = .01, 
I = .1, R = 200 nm, R 1 = 180 nm, R 2 = 220 nm, M = 8 x 10 5 A/m (permalloy), a = .01, 
and 7 = 1.7 x 10 1 T s . Deviation of low-field switching rate in (a) from dashed line signals 
non-Arrhenius behavior. 

magnetization density and large exchange constants. This combination occurs naturally in 
certain ferrimagnets. One example is MgOFe203, which has an exchange length a factor of 
five larger than that of permalloy. There are many examples of such materials that have 
been prepared as polycrystalline thin films, and thus are soft magnets (i.e., have no or very 
weak magnetocrystalline anisotropies) . Such materials might prove useful for experimental 
studies of the phenomena described in this paper. 



VI. DISCUSSION 



A theory of magnetization reversal in thin micromagnetic rings has been presented. Such 
systems are distinguished by their lack of edges or corners where nucleation is easily initiated, 
leading to greater stability of magnetization configurations and facilitating comparison of 
theory to experiment. 

By utilizing a scaling analysis [19] that uncovers a separation of energy scales in the thin 
film limit, we are able to retain leading-order terms that allow for an analytic solution of 
the relevant magnetization configurations in the thermally-induced reversal problem. The 
discarded terms, in particular those corresponding to nonlocal magnetostatic energy con- 
tributions, are shown to contribute no more than O(10%) to the energy over most values 
in the (£, h) phase plane. A complete solution that takes into account all terms must be 
numerical, and is planned for future work. 
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Nevertheless, an analytic solution is highly useful and can uncover information that may 
be difficult to extract from a numerical one. In particular, we predict an unusual transition 
from Arrhenius to non-Arrhenius activation behavior (cf. Fig. [7J). Our analysis suggests that 
such a transition should now be observable experimentally, by varying the externally applied 
magnetic field for rings of fixed size. A clear signature of such a transition would be the 
observation of a crossover from Arrhenius to non-Arrhenius behavior as field varies, as seen 
in Fig. [7| Because this requires measurement of the prefactor, such an observation would 
require numerous runs where reversal occurs. 

Arrhenius behavior of magnetic reversal has already been found in several systems and 
geometries. In Q], measurements of switching field and waiting times on nearly spherical Ni, 
Co, and Dy nanoparticles found an activation volume close to the particle volume, indicating 
a uniform magnetization reversal (analogous to the constant saddle case here) and confirming 
the Neel-Brown theory for these systems. In contrast, measurements on Ni wires with 
diameters 40-100 nm revealed an activation volume considerably smaller than the particle 
volume, indicating a nonuniform transition state (analogous to our instanton saddle). Here, 
too, Arrhenius switching behavior was found. But wouldn't the arguments given above imply 
that one should see non-Arrhenius behavior for these wires? No, because here the (roughly) 
cylindrical geometry with end caps (and consequently the relevant boundary conditions on 
the magnetization), lead to the absence of a uniform translation or rotational symmetry 
present, and therefore no zero mode leading to a temperature-dependent prefactor. In fact, 
these observations support the robustness of our conclusions for a wider variety of cases than 
considered in this paper. We will discuss this further below. 

Before doing that, however, we wish to suggest a second experimental test that may 
be easier to conduct: this is to measure the dependence of the activation energy on mean 
radius R for rings of identical composition. On the Arrhenius side of the transition, where 
magnetization reversal proceeds via a uniform rotation of the magnetization, the activation 
barrier scales linearly with the ring size. However, on the non-Arrhenius side, where the 
instanton state governs the reversal, the activation barrier is almost independent of ring size 
(see Fig. EJ). In this set of measurements, one may need to alter the applied field as ring size 
varies to keep the system on one or the other side of the transition, given that the critical 
field depends on R (cf. Eq. (jT2)). 

How robust are our predictions of a transition in activation behavior, and in particular, 
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can the neglected energy contributions wash out or obscure the transition? It is indeed 
possible, perhaps likely, that the details of the transition close to the critical field (or cir- 
cumference if field is fixed) are sensitive to these terms. In particular, the second-order 
nature of the transition, and the corresponding divergence of the prefactor (cf. Fig.|5J), could 
disappear. Inclusion of the magnetostatic terms could even in principle change the transi- 
tion from second- to first-order, with a jump replacing the divergence in the prefactor. Such 
first-order transitions have been predicted to occur in thermally-induced conductance jumps 
in monovalent metallic nanowires 

However, our central prediction, a transition from Arrhenius to non-Arrhenius activation 
behavior, should be robust because it is due to something much more fundamental: a ro- 
tationally invariant transition state (our 'constant' state <fi = cos _1 (— h)) at high fields and 
a rotationally non- invariant state (our instanton state ©) at low fields, with the crossover 
determined primarily through a competition between the shape anisotropy arising from 
magnetostatic forces and the Zeeman energy arising from the external field. In fact, the 
discussion in Sec. IV C 21 leads to the conclusion that the appearance at lower fields of any 
rotationally non-invariant state should give non-Arrhenius switching behavior in the ring 
geometry. Experimentally, what is then required is a symmetric enough ring so that the 
'domain wall' part of the transition state (centered at Sq in our instanton solution) has more 
or less equal probability of nucleating anyplace along the ring. This 'Goldstone mode', aris- 
ing from the rotationally invariant geometry, is ultimately where the non-Arrhenius factor 
comes from. Although the size constraints on the ring parameters leading to the specific 
instanton solution (JHJ) are difficult to realize at the present time, the generality of the basic 
physical features determining the transition should lead to the predicted crossover from Ar- 
rhenius to non-Arrhenius behavior in at least some ring geometries that are outside of the 
scaling regime considered here. 
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